BaNiF4: an electric field-switchable weak antiferromagnet 



o 
o 

(N 

cn 
(N 



X3 

o 

> 
O 
O 

O 

o 

-I— > 



I 

o 
o 



X 
S3 




Claude Ederer and Nicola A. Spaldin 
Materials Department, University of California, Santa Barbara, CA 93106-5050, U.S.A^ 

(Dated: February 6, 2008) 

We show that in the antiferromagnetic ferroelectric BaNiF4 the Dzyaloshinskii-Moriya interac- 
tion leads to a small canting of the magnetic moments, away from the easy axis, resulting in a 
noncoUinear magnetic structure. The canting corresponds to an additional "weak" antiferromag- 
netic order parameter whose orientation is determined by the polar structural distortion and can be 
reversed by switching the ferroelectric polarization with an electric field. Our results point the way 
to a more general coupling mechanism between structural distortions and magnetic order parame- 
ters in magnetoelectric multiferroics, which can be exploited in the design of electric field-switchable 
magnets. 



There is great interest in magnetic ferroelectrics, which 
primarily stems from the possible cross-correlations be- 
tween their magnetic and dielectric properties 0. Such 
magnetoelectric coupling facilitates the manipulation of 
magnetic properties using an electric field and vice versa. 
Various intriguing effects have been reported so far, rang- 
ing from a magnetic field dependence of the dielectric 
permittivity in various materials including YMnOa ^l, to 
magnetic phase control using an electric field in HoMnOs 
0, and the reorientation of electric polarization by a 
magnetic field in TbMnOs 0. 

Of particular appeal, both from the viewpoint of fun- 
damental science and also with respect to possible ap- 
plications in digital memory technologies, is the switch- 
ing of magnetic domains by an electric field and vice 
versa. Here, the field is used to push the system into 
a different realization of the same ground state phase. 
The resulting state continues to be stable even when the 
field is removed, thus exhibiting the basic requirement 
for nonvolatile data storage. Such magnetoelectric do- 
main switching has been demonstrated for multiferroic 
NiaByOial, where an electric field can be used to ro- 
tate the magnetization by ±90° and a magnetic field can 
switch the polarization into any possible domain config- 
uration j^. Nevertheless, despite the recent revival of 
interest in magnetoelectric phenomena, this has, to our 
knowledge, remained the only report of magnetoelectric 
domain switching in a single-phase multiferroic. In many 
multiferroics different mechanisms are driving the mag- 
netic and ferroelectric orders respectively, resulting in 
only weak coupling between the two order parameters, 
and disfavoring the possibility of magnetoelectric domain 
switching. However, electric field-induced magnetization 
switching has been observed recently in a nano-composite 
consisting of ferrimagnetic CoFe204 pillars embedded in 
an antiferromagnetic and ferroelectric BiFeOa matrix . 

In this letter we introduce a mechanism in which the 
magnetic order is induced by the ferroelectricity and is 
therefore intimately coupled to the polarization, allowing 
controlled switching of the magnetic order parameter by 
switching the electric polarization using an electric field. 
We demonstrate this approach using as an example the 



FIG. 1: (Color online) Projection of the BaMF4 structure 
along the a axis. The M cations are surrounded by fluorine 
octahedra, which form puckered sheets perpendicular to the 
h axis, separated by similar sheets of Ba cations. Note that 
adjacent sheets are shifted relative to each other by half a 
lattice constant along the a direction. 



magnetic ferroelectric BaNiF4, which is representative of 
a whole class of multiferroic barium fluorides. We show 
by using both first-principles calculations and symmetry 
analysis that the polar distortion in this material gives 
rise to an additional "weak" antiferromagnetic order pa- 
rameter which is reversed when the ferroelectric polar- 
ization is reversed. We suggest the possible experimental 
observation of this electric field-switchable weak antifer- 
romagnetism and discuss the general implications of our 
results for magnetoelectric switching phenomena. 

BaNiF4 is representative of the isostructural family 
of barium fluorides with the chemical formula BaAfF4, 
where M can be Mn, Fe, Co, Ni, Zn, or Mg. These com- 
pounds were first synthesized in 1968 0,13' crys- 
tallize in a base-centered orthorhombic structure with 
space group Cmc2i Q which is shown in Fig. ^ 

BaMnF4 undergoes an additional structural phase tran- 
sition at around 250 K (see Ref. 0). Ferroelectric switch- 
ing has been reported for M=Co, Ni, Zn, and Mg but 
not for M=Mn and Fe [13. The systems with M=Mn, 
Fe, Co and Ni order antiferromagnetically at Neel tem- 
peratures around 20-70 K. The magnetic structure de- 
rived experimentally for the systems with M=Mn, Fe, 
Ni is shown in Fig. |2K The magnetic unit cell is 

doubled compared to the chemical unit cell, and con- 
tains four magnetic M cations, which are arranged in 




FIG. 2: (Color online) Magnetic structures (not to scale) of 
BaNiF4 derived from the experimental observations (a) and 
from our calculation including spin-orbit coupling (b). Gray 
lines outline the conventional orthorhombic unit cell, green 
lines show the puckered sheets perpendicular to the b direc- 
tion, numbers in (a) indicate the 4 magnetic ions in the unit 
ceU. 



sheets perpendicular to the b axis. Within each sheet 
the cations form a puckered rectangular grid, with the 
magnetic moments of neighboring cations oriented an- 
tiferromagnetically, and all moments aligned parallel to 
the b axis (c axis for M=Co 12]). The coupling between 
different sheets is weak, leading to low magnetic ordering 
temperatures and pronounced two-dimensional behavior 
[l^. The corresponding magnetic space group is Pa^i 
(Pa2i for BaCoF4). 

In this letter we focus on the intriguing magnetoelectric 
and magnetostructural properties of BaNiF4. A detailed 
comparative study of the structural, electronic, and mag- 
netic properties of the whole series of multiferroic com- 
pounds will be reported elsewhere All calculations 
presented in this work are performed using the Vienna 
Ab-Initio Simulation Package (VASP) ^14\ employ- 
ing the projector-augmented wave method 15, 16J. To 
account for the strong Coulomb interaction between the 
localized d electrons of the transition metal cations we 
use the LSDA+C/ method 0. For the Hubbard U and 
intra-atomic exchange parameter J we use typical val- 
ues of C/=4eV and J=leV, and we check the sensitiv- 
ity of our results with respect to variations of both U 
and J. We use a plane- wave energy cutoff of 450 eV 
and a F-centered 2x4x2 k-point mesh (divisions with 
respect to the monoclinic basis vectors of the magnetic 
unit cell), which is sufficient to obtain converged results 
for all quantities under consideration. 

First we perform a full structural optimization of 
BaNiF4 within the experimentally observed Cmc2i sym- 
metry; the structural parameters we obtain agree well 
with available experimental data |l3l| . To confirm the 
experimentally reported magnetic structure we then cal- 
culate the energy differences for different relative orien- 
tations of the four magnetic sublattices corresponding to 
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FIG. 3: (Color online) Dependence of canting angle on the 
LSDA-I-C/ parameters. Black circles correspond to J=0 and 
red squares to J—1 eV. 

the four magnetic ions in the unit cell and extract the 
Heisenberg nearest neighbor coupling constants. In these 
calculations the spin-orbit coupling is neglected. We ob- 
tain strong antiferromagnetic nearest neighbor coupling 
within the buckled planes, Ja=4.4meV, Jc=3.3meV |33, 
consistent with the experimentally observed structure. 
However, when spin-orbit coupling is included in the cal- 
culation, we observe that the collinear spin-configuration 
shown in Fig. Et, with all spins aligned along the b di- 
rection, is unstable, and that instead the magnetic mo- 
ments assume a noncoUinear configuration where all spins 
are slightly tilted towards the ±c direction as shown in 
Fig. [2t>- The tilting angle is about 3°. 

Figure 121 shows the U and J dependence of the canting 
angle between the magnetic moments of the Ni cations 
and the b direction. In order to obtain local magnetic 
moments we integrate the spin-density within spheres 
of radius 1.2 A centered at the Ni sites. Our test cal- 
culations show that for this radius the integrated spin- 
density adopts a saturated value of about 1.7 hb- It can 
be seen that within the physically reasonable range of 
U—2-6eV there is only a moderate U dependence of the 
canting angle, whereas it is slightly suppressed by de- 
creasing J. This J dependence is a consequence of the 
spin-non-diagonal elements introduced into the effective 
potential by the LSDA+J7 energy in its most general form 
[13], which is rotationally invariant with respect to both 
orbital and spin quantum numbers. It follows from Fig. El 
that within the physically reasonable parameter range of 
C/=2-6eV and J=0-leV the canting angle is 2-3°. 

The small tilting of the magnetic ions can be explained 
by the antisymmetric exchange or Dzyaloshinskii-Moriya 
(DM) interaction [l^], = dij-{siXSj), where s, is the 
spin of ion i and dij is the coupling vector correspond- 
ing to the antisymmetric exchange interaction between 
ions i and j. This interaction occurs only for certain low 
symmetries and can give rise to "weak ferromagnetism" 
in otherwise antiferromagnetic materials, where a canting 
of the two (or more) sublattice magnetizations away from 
the ideal collinear antiferromagnetic orientation gives rise 
to a small net magnetization 0,113. In the case of 
BaNiF4 the magnetic space group does not allow the oc- 
currence of weak ferromagnetism but nevertheless there 
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is a nonzero DM interaction dc between magnetic near- 
est neighbors along the c direction [s^l , whereas the DM 
interaction between neighboring cations along the a di- 
rection vanishes by symmetry. 

The canting due to the DM interaction leads to a 
"weak" antiferromagnetic order parameter — S1 + S2 — 
S3 — S4, in addition to the experimentally observed (pri- 
mary) antiferromagnetic order parameter Lab = si — S2 — 
S3 + S4. On a macroscopic level the DM interaction leads 
to a coupling between Lab and L^ of the form: 

EY^,,,^D-{LabxL,) , (1) 

where D = dc/2. If we neglect the magnetic single-ion 
anisotropy contribution to the total energy, the canting 
angle is given by a « D/AJc, from which we obtain a 
value for D of about 0.7 meV. 

Note that a term of the form can also be caused 
by different orientations of the easy axes for the mag- 
netic moments on different sites. This mechanism is also 
allowed within the crystallographic and magnetic sym- 
metry of BaNiF4, in addition to the DM interaction de- 
scribed above. Since the following analysis of the magne- 
toelectric domain switching is independent of the actual 
microscopic mechanism leading to the term we do 
not separate these two effects further. 

In order to analyze the possibility of magnetoelectric 
domain switching in BaNiF4 we need to know the sym- 
metry of the corresponding "prototype" phase [i^, 0| ■ 
Although the paraelectric states of the BaMF4 systems 
are not accessible experimentally, since all crystals melt 
before undergoing a ferroelectric phase transition, a non- 
polar reference structure has been suggested and used to 
discuss the ferroelectric switching properties of these sys- 
tems ll^. This structure has the nonpolar space group 
Cmcm and can be obtained from the ground state Cmc2i 
structure by enforcing a mirror symmetry perpendicu- 
lar to the c axis. Within this symmetry no canting of 
the magnetic moments is allowed and the resulting mag- 
netic order corresponds to the coUinear spin arrangement 
shown in Fig.|21i with the magnetic space group Pa^i/m! . 
The absence of canting in the centrosymmetric reference 
structure indicates an intimate connection between the 
structural distortion and the weak antiferromagnetic or- 
der parameter L^. 

In BaNiF4 only two different structural domains are 
possible, corresponding to the two opposite orientations 
of the ferroelectric polarization along ±c. Indeed, if we 
invert the polar mode in our calculation, we observe that 
for fixed orientation of Lab the orientation of the weak 
antiferromagnetic order parameter is determined by the 
orientation of the polarization P (see Fig.0J. This means 
that a reversal of P by an electric field is accompanied by 
a reversal of Lc, if during the polarization switching the 
primary antiferromagnetic order parameter Lab is pre- 
served. Energetically, a reversal of Lc is much more fa- 
vorable than reversal of Lab, since the latter requires the 




FIG. 4: (Color online) Magnetoelectric switching in BaNiF4. 
The weak antiferromagnetic order parameter Lc is coupled to 
the polarization P. Reversal of the polarization from (a) to 
(b) leads to a reversal of the canting of the magnetic moments 
and thus to a reversal of Lc- 



magnetic moments to rotate by ^ 180 ° through the hard 
magnetic axis, whereas the reversal of L^ requires only 
a slight reorientation of the magnetic moments through 
the easy axis. We therefore propose that the ferroelectric 
switching in BaNiF4 will be accompanied by a reversal 
of the weak antiferromagnetic order parameter. 

A canted antiferromagnetic structure caused by the an- 
tisymmetric exchange interaction has also been proposed 
in the seminal paper by Moriya 0^ for CuCl2 -21120, and 
has subsequently been confirmed by several experimental 
methods including neutron scattering 23] , antiferromag- 
netic resonance measurements j2Jj, and NMR [2^. This 
shows that a "hidden" order parameter such as Lc can 
be detected experimentally. However, in order to observe 
the predicted reversal of Lc by an electric field, an exper- 
imental technique is required that allows the distinction 
between different antiferromagnetic 180° domains. One 
such technique, which in addition allows imaging of the 
antiferromagnetic domain topology, is second-harmonic 
generation (SHG) [2^, which, in fact, has already suc- 
cessfully detected weak antiferromagnetism in YMnOa 
[27I 1. We therefore hope that our work will stimulate ex- 
perimental efforts to measure the predicted magnetoelec- 
tric domain switching in BaNiF4 and related materials. 

Magnetic order that is induced by ferroelectricity via 
the magnetoelectric effect has been suggested for the re- 
lated system BaMnF4, based on macroscopic symme- 
try considerations [2^. In BaMnF4 an additional struc- 
tural phase transition lowers the symmetry so as to al- 
low the linear magnetoelectric effect, which is symmetry- 
forbidden in BaNiF4. Nevertheless, the microscopic 
mechanism governing magnetoelectric behavior could be 
the same in both systems (although a rigorous micro- 
scopic study is required to show this unambiguously). 
The weak antiferromagnetism in BaNiF4 would then be 
classified as "anti-magnetoelectric effect" mediated 
by the DM interaction. A related scenario has also been 
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reported recently for multiferroic BiFeOa , where the DM 
interaction leads to weak ferromagnetism . Although 
the linear magnetoelectric effect is symmetry-allowed in 
this system |34| . first principles calculations show that 
the weak ferromagnetic moment in BiFcOs is in fact not 
induced by the ferroelectricity, but is coupled to a non- 
polar structural mode [2oj . which exemplifies that sym- 
metry considerations alone cannot provide full insight 
into the exact coupling mechanism of magnetoelectric 
systems. 

Since the DM interaction can only be nonzero if the 
midpoint between two magnetic ions is not an inversion 
center, the occurrence of "weak" magnetic order due to 
spin-canting should be a rather common phenomenon in 
magnetic ferroelectrics. The resulting order parameter 
can in general be ferro- or antiferromagnetic and can be 
reversed by reversing the structural mode that is respon- 
sible for the nonzero DM interaction. If this mode is 
polar, it can be reversed by an electric field, thus leading 
to an electric field-switchable magnetic order parameter. 
Our results for BaNiF4 and our earlier study of BiFcOa 
show that the amplitude of the corresponding order pa- 
rameters are typically of the order of 0.1 /ie per magnetic 
ion, which is about 10 % of a typical "strong" magnetic 
order parameter (~ 1 /xb per magnetic ion) , which is easy 
to measure with modern experimental techniques. 

Finally, we compare this with a somewhat converse 
effect, which is currently attracting a lot of attention, 
namely the appearance of a small electric polarization 
in materials where the magnetic order breaks the inver- 
sion symmetry of the system. The resulting polariza- 
tion is then intimately coupled to the magnetic properties 
0, l29j | . Although the underlying microscopic mechanism 
is not yet understood, the important role that symme- 
try plays for this effect is similar to the case of the weak 
magnetic order discussed in this letter. The magnitude of 
the resulting electric polarization is around 0.01 /iC/cm^ 
0,113 , which is two to three orders of magnitudes smaller 
than the spontaneous polarization in typical "strong" fer- 
roelectrics (~1-10 /iC/cm^). From this it becomes clear 
that the "weak" magnetic order that is caused by struc- 
tural distortions is in comparison a rather large effect 
and therefore provides a very promising route for achiev- 
ing practical magnetoelectric coupling effects. In addi- 
tion, since the weak magnetic order parameter and the 
ferroelectric polarization are coupled by symmetry, it is 
not necessary that the primary magnetic and ferroelec- 
tric order temperatures are close together. The weak 
magnetic order appears at the Neel temperature but is 
in fact caused by the structural distortions, which can 
"freeze in" at much higher temperatures. 

In summary, we have shown that the multiferroic sys- 
tem BaNiF4 exhibits an additional "weak" antiferromag- 
netic order parameter that is coupled to the polar dis- 
tortion, and can be reversed by reversing the ferroelec- 
tric polarization using an electric field. We propose this 



mechanism as a promising route for achieving magneto- 
electric switching behavior in ferroelectric magnets. 
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